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Using cooling for SU{2) lattice configurations, purely Abelian constant magnetic field configu- 
rations were left over after the annihilation of constituents that formed metastable Q = con- 
figurations. These so-called Dirac sheet configurations were found to be stable if emerging from 
the confined phase, close to the deconfinement phase transition, provided their Polyakov loop was 
sufficiently non-trivial. Ifere we show how this is related to the notion of marginal stability of the 
appropriate constant magnetic field configurations. We find a perfect agreement between the ana- 
lytic prediction for the dependence of stability on the value of the Polyakov loop (the holonomy) in 
a finite volume and the numerical results studied on a finite lattice in the context of the Dirac sheet 
configurations. 
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I. INTRODUCTION 

In lattice gauge theory cooling is used to remove the 
high frequency fluctuations to be left with classical so- 
lutions 0, 12I • This allows one to extract the underlying 
topological content of the gauge field configurations and 
determine to what extent instantons have a role to play. 
It is known that when using the ordinary Wilson action, 
the lattice artefacts are such that one can further lower 
the action by reducing the size of the instantons (whereas 
in the continuum the classical action does not depend on 
the size) . Ultimately the instanton falls through the lat- 
tice and in general one relaxes to the trivial minimum 
with zero action. However, at finite temperature, when 
the Polyakov loop away from the instanton is non-trivial, 
the relevant instanton (called a caloron) actually consists 
of n constituents for SU(n) These can be shown to 

be 't Hooft-Polyakov (BPS) monopoles when identifying 
Ao with the (adjoint) Higgs field. From the Euclidean 
four-dimcnsional point of view, due to the self-duality 
of the gauge field, these are dyons with their magnetic 
charge equal to their electric charge, with overall electric 
and magnetic neutrality. 

Under cooling in the confined phase, due to the dis- 
creteness artefacts of the Wilson action, these con- 
stituents will attract and approach each other. When 
they are no longer visible as separate entities, the solu- 
tions behave like ordinary instantons localized in space 
and time. The distance between the constituents is (for 
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SU(2)) given by irp^ /b, where b is the inverse temperature 
(the period in the euclidean time direction). Another 
possibility is the annihilation of dyons and antidyons left 
over from different caloron and anticaloron solutions. As 
a result, with an action near the one-instanton action, 
a metastable configuration can be either a dyon-dyon 
pair that shrinks and falls through the lattice or a dyon- 
antidyon pair that finally annihilates |3| • Sometimes this 
annihilation process leaves behind a constant Abelian 
magnetic field, which subsequently turns out to be stable 
or unstable under further cooling |^ , strongly correlated 
to the asymptotic value of the Polyakov loop (the holon- 
omy) which has been acquired in this stage of cooling. In 
the deconfined phase no dyonic structure was observable 
under cooling. The Polyakov loop remains always close 
to its trivial value but quasi-constant magnetic field con- 
figurations were seen to emerge as well, although they 
never happened to be stable. In this paper we present an 
explanation for these observations. 

II. CONSTANT MAGNETIC FIELDS 

It is well-know that Abelian constant magnetic fields 
are embedded solutions of the (non- Abelian) equations 
of motion. They tend to be unstable, due to the self- 
coupling of the gauge fields 7], which formed the basis 
for the studies of the so-called Copenhagen vacuum pic- 
ture 

In the four-dimcnsional context a constant field is sta- 
ble if it is self-dual 0, ^3 ■ Here we will be interested in 
the degenerate case with magnetic, but no electric flux 
and periodic boundary conditions (the general case al- 
lows for center fluxes, but requires twisted boundary con- 
ditions ITU). For SU{2) these gauge fields are Abelian 
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and there is the freedom of adding to it a constant 
Abehan vector potential, which does not change the field 
strength, F^^, = 7rir3n^^/(L^Li,). This field strength 
is unique up to a constant gauge rotation, and n^i, is 
an integer (even in the case of periodic boundary condi- 
tions) antisymmetric tensor, fixed by flux quantization. 
In the degenerate case n^j^ has two non-zero eigenvec- 
tors, and computing the gauge-invariant Polyakov-loop 
observables in this subspace it is easily seen that no trans- 
lation invariance holds. Adding a constant Abelian vec- 
tor potential can consequently be absorbed by a trans- 
lation and can therefore not affect the fluctuation spec- 
trum. But in this degenerate case there are also two zero 
eigenvectors, and the vector potential is invariant under 
translations in this subspace. Its Polyakov loops label 
the gauge invariant parameters on which the fluctuation 
spectrum do depend! 

It had been found 0, 0| that on a symmetric torus 
there was one class of constant magnetic field solutions 
that for a certain range of values of the Polyakov loop 
were stable. This example, involving the smallest possi- 
ble non-zero magnetic field, possesses non-trivial center 
flux and requires twisted boundary conditions. Therefore 
it could not explain the flndings of Ref. (Qj. However, at 
finite temperature involving a non-symmetric box, more 
room exists to obtain stable constant magnetic field so- 
lutions. 

III. THE FLUCTUATION SPECTRUM 

For SU{2) all constant curvature solutions in a finite 
box have been classified. Also the spectrum of fluctua- 
tions has been calculated For the "charged" isospin 
components, in the subspace of non-zero eigenvectors of 
n^iy, the problem is equivalent to that of Landau levels. 
The eigenf unctions are described by O functions to in- 
corporate the boundary conditions. In the subspace of 
zero eigenvectors one simply has plane waves, with prop- 
erly discretized momenta. These momenta are, however, 
shifted due to the constant vector potential which deter- 
mines the Polyakov loops in this subspace, thereby obvi- 
ously modifying the fluctuation spectrum. The "neutral" 
isospin component is described by ordinary plane waves. 

The following gauge field for SU{2) gives the most gen- 
eral solution with constant field strength on a torus [iJI 

-^u{x) = ^i{-Tmi,f_,Xfj,/Lfj, + Cu)t3/Li, . (1) 
It is periodic up to the gauge transformation 

A,{x + fiL^) ^ n^{x){AAx) + d,)n-'{x) (2) 

where ft is the unit vector in the fi direction and 

^f^ix) = exp{^iTTXi,nt,f,T3/L„) . (3) 

With 71^1^ even, these Abelian boundary conditions are, 
however, gauge equivalent (in general by a non- Abelian 
gauge transformation) to periodic boundary conditions 



(as long as Q = 0). Following Ref. we assume Lq = 
Lt = b, Li — L2 = — Ls- The data can in all cases 
be interpreted in terms of a (nearly) constant magnetic 
field with noi, — —n^Q — and rh ~ (0,0,2), where 
rrii — ^SijkUjk- Therefore we compute the fiuctuation 
eigenvalues for this case (compare Ref. [11 El El) 

A± = 47r(2n+ 1±2)/L2 

+ {2nq + Co)y^, (4) 
Ao = (2^fc^/L^)2 . 

The multiplicities are 4 for A± and 2 for Aq, with all 
quantum numbers {n,p, q, k^) integer (but n > 0). 

As argued above the spectrum depends on the constant 
Abelian gauge field described by the constants Co and C3 . 
These are only defined modulo 27r, as a shift over 27r is 
related to a gauge transformation that shifts the relevant 
momenta by one unit. The Polyakov-loop observables 
are given by 

= iTrexp(iC^T3/2) = cos(C^/2) , = 0, 3 . (5) 

Note that these are anti-periodic under a shift over 27r, 
whereas the fluctuation spectrum is periodic. This is sim- 
ply because the fluctuations involve fields in the adjoint 
representation, whereas the Polyakov loop is in the fun- 
damental representation. Indeed Pq and P| , relevant for 
the adjoint representation, are periodic under a shift of 
Co and C3 over 27r. 

From the lattice data it is clear that P3 = 1, and we 
can put C3 — 0, as well 'as n — p — q — Q (we may 
restrict |Co| < tt) to find the lowest eigenvalue A_ = 
—^■k/LI + Cq/LI to be negative unless the Polyakov loop 
is sufficiently non-trivial (Po = ±1 being associated to a 
trivial Polyakov loop). The lowest eigenvalue is positive 
when 

Lt/Ls < 

(6) 

|Pq| = cos(Co/2) < cos{^/^Lt/Ls) . 

We see that these conditions cannot be satisfied when 
Lt = Lg and the finite temperature situation [b = Lt < 
Lg) is essential for providing the opportunity of stability. 

This stability was called marginal, because one can 
change Co without changing the classical action. Thus 
nothing prevents us to bring Co close to 0, where the 
lowest eigenvalue A_ turns negative. Under the cooling 
there is no reason for Co to change, as one can easily 
show that the degeneracy of the action as a function of Cq 
survives on the lattice. This then explains the stability 
of these constant gauge field configurations, provided the 
two conditions of Eq. ijHJl are satisfied. 

IV. COMPARISON WITH LATTICE DATA 

In Ref. SU{2) gauge theory in four-dimensional Eu- 
clidean space was considered on an asymmetric lattice 
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with periodic boundary conditions in all four directions. 
The respective ensembles of configurations have been cre- 
ated by heat-bath Monte Carlo using the standard Wil- 
son plaquette action. The lattice size was x Nt with 
the spatial extension Ng = 8, 10, 12, 16, 20 and with tem- 
poral extension Nt = 4, i.e. b = 4a and Lg = aNs with 
a the lattice spacing) . For Nt = A the model is known to 
undergo the deconfineme nt p hase transition at the crit- 
ical coupling Pc — 2.299 0|. In Ref. two ensembles 
with (3i = 2.2 < f3c and /32 = 2.4 > (5^ were generated. 

The equilibrium field configurations in both ensembles 
have been cooled by iterative minimization of the Wilson 
action with the focus at the structure of selfdual caloron 
solutions. In addition, Dirac sheet (DS) events were ob- 
served at the very last stages of cooling, applying a stop- 
ping criterion which selects action plateaux in the inter- 
val S < 0.6 ^inst- In confinement, approxmately 7 % (at 
Ns = 8,10,12), 5 % {Ns = 16) and 3 % {Ns = 20) of 
equilibrium configurations have turned into these purely 
magnetic configurations, whereas in the deconfinement 
phase the yield was 5 • • • 18 % 0. The action values 
were found close to {Nt/Ns)Sj^ characteristic for con- 
stant Abelian magnetic flux [ig of size Att periodically 
closed along one of the spatial directions. Although the 
action showed the same dependence on the lattice exten- 
sions Nt and Ng, supporting the common interpretation 
as (almost) homogeneous magnetic flux, the configura- 
tions were unstable when derived from the deconfined 
phase and partly stable in the case of the confinement 
phase. In the case of confinement, the issue of stability 
vs. instability was strongly correlated to the value of the 
temporal Polyakov line (holonomy) Pq. This is shown 
in Figure ^ It presents a set of scatterplots (each for 
another spatial size Ng) where each DS event is charac- 
terized by two entries: (sm™, Pq) and {smax, Po)- The 
values Smin and Smax express the action density at sites 
where it is minimal and maximal, respectively. If these 
values differ, the configuration is bound to decay to the 
trivial vacuum. Provided the holonomy remains suffi- 
ciently far from trivial, we find only DS events which 
consist of a highly homogeneous Abelian magnetic flux 
signaled by Smin — Smax- This case is tantamont to ab- 



solute stability under further cooling. In contrast to this, 
when the holonomy was close to the trivial one (P « ±1) 
the Abelian magnetic fluxes happened not to be homo- 
geneous {srnini 7^ Smax) and provcd to be unstable un- 
der further cooling. The critical value of the holonomy, 
I Pol = cos{^ Lt / Lg), limiting the region of stability as 
given by the second condition in Eq. lO, is marked in 
Figure^by horizontal lines. No deviations from the pre- 
dicted (in)stability are seen. 

V. CONCLUSIONS 

Purely Abelian constant magnetic field configurations 
were observed |^, randomly emerging from the process 
of cooling-down equilibrium lattice fields representing the 
confined and deconfined phases of SU{2) gluodynamics. 
In the confined phase they were found to be absolutely 
stable provided their Polyakov loop was sufficiently non- 
trivial. We have shown here that this fact is related to 
the notion of marginal stability of the appropriate con- 
stant magnetic field configurations. We have found per- 
fect agreement between the analytically predicted depen- 
dence of stability on the value of the Polyakov loop (the 
holonomy) for the set of spatial lattice sizes that were 
studied in Ref. 0| and the numerical observations made 
there, separating stable from unstable Dirac sheet config- 
urations. The dependence on the geometry of the effect 
we found makes us believe we are dealing with a finite 
volume artefact. Nevertheless, it demonstrates that the 
influence of a background constant Aq (as manifested by 
the Polyakov loop) on the dynamics of the gauge field 
should not be ignored. The physical significance lies in 
the fact that the Polyakov loop is the order parameter for 
the confinement /deconflnement phase transition. A sim- 
ilar conclusion can be drawn from the caloron solutions 
with non-trivial holonomy. 
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FIG. 1: Correlation in the confined phase of 5(7(2) LGT between the holonomy Pq and the (in)stability of Dirac sheets which 
is indicated through the remaining inhomogeneity of the action density, Smax 5^ Smin (represented by circles and triangles, 
respectively). The limiting values of the holonomy, |-Po| = cos(-y7r^), beyond which constant magnetic fluxes become unstable 
are indicated by horizontal lines. The temporal size of the lattices is Nt = 4, and the DS events are shown for different spatial 
lattice size Ns = 8, 10, 12, 16, 20. 
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